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We study the insertion and value of metric observables in the (discrete) path integral formulation 
... of the Ponzano-Regge spinfoam model for 3d quantum gravity. In particular, we discuss the length 

OO ' spectrum and the relation between insertion of such B-observables and gauge fixing in the path 

' integral. 
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I. INTRODUCTION 



As it stands, it is technically extremely difficult to extract physical information from gravitational theories, 
both classical and quantum level. In both regimes, this process requires the definition of observables, which 
O should be both invariant under the symmetries of the theory, and nonetheless, physically relevant. Such 

observables are hard to come by, let alone compute. At the classical level, the majority of work has been 
^ [ done within the canonical setting. Only recently, has the road been opened to fully calculate observables 

within this formalism [l| , using the technique of partial observables . 

At the quantum level, however, almost nothing is known of physical geometric observables. Some progress 
has been made within both the canonical quantisation process and the path integral approach. In loop 
^ ' quantum gravity the geometric observables lie mostly in the domain of the kinematical Hilbert space 

, where the volume and area operators are defined and are diagonal on the space of states [3| • In the realm of 

covariant approaches to quantum gravity, a systematic analysis of the space of physical geometric observables 
has remained untouched. It is to redress this balance that we apply our resources hereafter. 

On the other hand, much progress has been made recently, within both the canonical and covariant 
OO ' formalisms, to evaluate matter observables [1, 0, 0, Q (both at the classical and quantum level). Within 

. the discrete path integral approach known as spin foams [9, 10], the effective matter observable have been 

shown to arise as the Feynman diagrams of a non-commutative quantum field theory [ll|, [l^] ■ Interestingly, 
matter has also been dealt with in the regime of group field theory [isj ; a theory which generates spin foams 
as Feynman diagrams. Recently, the non-commutative field theory encountered above has been shown to 
arise as a particular perturbation around a classical solution of the group field theory [l4|. 
; ^ ■ We shall also follow the route of spin foams to define our path integral. These are sums over configurations 

I with support on piecewise linear manifolds. The amplitudes are purely combinatorial and depend on the 

representation theory of the appropriate gauge group. They may be recast in a form which resembles the path 
integral of discretised (constrained) BF theories. It is in this form, that we will perform our calculations. 
In the case of 3d gravity, there are no constraints, and we are left with the simple BF term. Therefore, the 
possible observables are gauge invariant functions of the B variables and the holomonies. 

We shall pay most attention here to geometric variables such as the length squared, and the dihedral angle 
between edges. To get meaningful results, one must gauge-fix the amplitude and we pay special attention 
here to this process. In the usual scenario [lH|, gauge fixing is the process of picking maximal tree with the 
relevant discrete structures and trivialising the variables on those edges of the trees. In general, a Fadeev- 
Popov determinant will arise, which we must take into account. We invoke this process here, also. To extract 
more information from the observables, we introduce normalised "states" which serve to pick out the edge 
of the triangulation upon which we intend to evaluate the observable. 
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Another motivation, is that B functionals, occur frequently, when we rewrite theories, such as Yang-MiUs 
and 4d quantum gravity, in terms of constrained BF theories [lol. [l6j. Thus, understanding their nature in 
this simphfied setting will pave the way to dealing with them more generally. 

The outline of the paper is a follows. In Section |lTl we outline the passage from the classical continuum 
theory, to the discrete path integral approach. In succeeding subsections, we recover the SU(2) Ponzano- 
Regge model from the S0(3) theory, and discuss the ambiguity in the discretisation of the B field. We 
then proceed to evaluate the expectation value of B observables in Section IIIII We reserve a more thorough 
investigation of the gauge-fixing procedure for Section IIVI and subsequently present our conclusions and 
outlook. 



II. DISCRETE PATH INTEGRAL FOR THE PONZANO-REGGE MODEL 



Before we commence to describe the main import of our work, we shall lay out succinctly the particular 
quantum gravity theory which provides the framework within which we shall calculate. For our purposes, 
it is sufficient to restrict ourselves to a 3d Riemannian theory of gravity without cosmological constant. To 
pass to a quantum theory, we shall write the action in \st order form as an 5u(2) gauge theory: 

S[B,uj] = ( tr(BAFM). (1) 

JM 

is a closed manifold. The triad _B is a su(2)-valued 1-form, from which the metric is reconstructed as 
g^i, = B^^Blrjij. The parallel transport on the manifold is given by the su(2)-valued connection uj and its 
curvature F[uj] — dcu + lu A oj is an su(2)-valued 2- form, tr is the trace over the Lie algebra. The classical 
equations of motion are: 

d^B = dB + [lu, B] = 0, F[w] = (2) 

where [•, •] is the Lie bracket on the algebra. These equations impose that the connection is torsion-free and 
flat. Solving for the connection in terms of the triad returns us to the standard 2nd order metric theory for 
3d Riemannian gravity. ?,j are su(2) indices and rjij is the flat metric on su(2), while ^, are 'space-time' 
indices. 

This action has several symmetries; namely, it is invariant under rotations, translations and diffeomor- 
phisms.^ These symmetries require proper gauge- fixing. Nevertheless, we can formally define the partition 
function"^ : 

Zm^ I VBVuj e'^^^^'^l (6) 



^ At the continuum level, we write down the gauge symmetries of this action: 



Rotation ^ ^k~^Bk~ ^ parametrised by A; € SU{2), (3) 



Translation s q d i j j, parametrised by </< G su(2). (4) 

The translation symmetry holds due to the Bianchi identity, d^F[ui\ = 0, provided </< = on the boundaries dM (i.e. the 
translation symmetry does not extend to the boundary). Remarkably, infinitesimal isu(2) symmetry is equivalent on-shell to 
difli'eomorphism symmetry provided det(i3) ^ 0. We see this by the following table of transformations: 

S^uj = d^p S^B = [B,p], 

5> = 0, 5jB = d^0, (5) 

5f^ = ^ + ^ + (FM), 5f B = 5«„ B + Sf^^ B + t^{d^B), 

where the diffeomorphism, parametrised by a vector acts like a Lie derivative: 5^ = di^ + %^d.^ The symmetry becomes 
manifest on-shell if we choose p = i^ui and <f> = i^B. 
^ This should not be confused with 3d Euclidean quantum gravity where the coefficient of the action in the partition function 
is —1 rather than i and may be regarded as an analytic continuation of the Lorentzian theory. 
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The triad field B is interpreted as a Lagrange multiplier enforcing the flatness of the connection: Z = 
J 'Duj6{F[uj])). As the BF theory action is topological, we may apply methods developed in topological 
quantum field theory to evaluate this partition function and appropriate gauge-invariant observables. This 
is not our aim here; we wish to study the link with discrete models of quantum gravity and compute 
the spectrum of geometrical operators (observables depending on the i3-field) in a discrete path integral 
formalism. 

With this target in mind, we replace the manifold by a simplicial manifold A of the same topology. 
Heuristically speaking, since the BF theory is topological and does not have any local degree of freedom, we 
do not expect to lose information in this replacement. In three space-time dimensions, we can triangulate 
any manifold. We label the 0-,l-,2- and 3-subsimplices as v,e,f and t, respectively. Another important 
constituent is the topological dual A* to the simplicial complex. We label sub-elements of this structure as 
v*,e*, f* and t* , respectively. The fields B and w are replaced by configurations which are distributional 
with support on subsimplices of A and its topological dual A*. We integrate these fields over the appropriate 
subsimplices. The definition of the integrated fields is: 



B 


B, 


= LB 


e SU(2), 




9e* 




e su(2), 


m - 




= Ue-edf 9e* 


e SU(2). 



The flatness constraint F[uj] = translates into the triviality of holonomies Ge = I (around closed loops, for 
trivial homotopy). At the discrete level, we will therefore replace the (5(F[a;]) constraint by d{Ge) constraints, 
with the discrete Bf, variables still playing the role of the Lagrange multipliers. Following [3], the action on 
the simplicial manifold then reads: 

eGA 

where the precise definition of the trace will be given later. Comparing this discrete action with the continuous 
BF action, we have simply replaced the curvature F by the holonomy G. The first remark concerns the 
definition of the holonomy Ge as the oriented product of the group elements ge* for all dual edges e* around 
the dual face (or plaquette) /* corresponding to the considered edge e. In order to define such a holonomy, 
one needs to choose a starting point along the loop. Changing the starting point amounts to acting on Ge by 
conjugation by some SU(2) group element k. This is compensated by rotating the _B-variable, Be — > kBek~^ . 
Therefore, as long as we restrict ourselves to considering gauge invariant observables of the B variables, such 
as the norm of Be, there is no issue. The second remark is about the difference between F[uj\ which lives in 
the Lie algebra su(2) and the holonomy Ge defined as a group clement in SU(2). Actually the original path 
integral derivation of the Ponzano-Regge model introduced the Lie algebra element Ze = log(Ge) G su(2) 
and considered the action iv{BeZe) Nevertheless, one faces the issue of the non-continuity of the log 
map (and the choice of a particular branch and so on). Therefore, we prefer to work with the formulation 
presented in Q which uses directly the group element Ge- 

The continuous gauge symmetries are broken, but there are residual discrete gauge symmetries: the 
rotation symmetry acts at the vertices v* E A* while the translation symmetry acts ai v £ A.* For details 



The gauge symmetries are: 



Rotation 



Translation 



9e* 
Be 

Qe* 
Be 



t(e*) 



<=*) 



Bek. 



v*(e) 



* ge 

Be 



■ U'J'Utie) 



■U!^''Usie) + l^°J"\<t>s(e)] 



parametrised by k^r € SU{2), (9) 



parametrised by <j>y G su(2), (10) 



where s(e*), t{e*) are the source and target vertices of e* , and v*{e) is the vertex where the holonomy around the face /* 
begins. It does not matter which vertex we choose as we may transfer to any vertex v* C df* by conjugation. Finally, 
= J^uj starting from v, in the limit where the edgelength is small , and C/J is a precise function of the holonomies given 
in Appendix iBl 
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on the action of these discrete gauge symmetries and their gauge fixing, we refer the interested reader to 0] ■ 
Finally, the discrete path integral reads: 

2A=n / d'B, n / d5e. e'^[^-^-l, (11) 

where the Lebesgue measure (PB^ on su(2) ^ and the Haar measure dge- on SU(2) are the natural choice 
for the discretised path integral measure. Due to the gauge invariance of the theory, to arrive at sensible 
results we must fix the gauge, [7]. To do so, we pick maximal trees of edges T G A and T* € A*.^ One then 
fixes Be = for all e € T and ^e* = I for all e* G T*. Mercifully, the Fadeev-Popov determinant for this 
gauge evaluates to unity. We will discuss this gauge fixing in more detail in the Section IIVI 



A. Recovering the Ponzano-Regge model 



Our first task is to evaluate the discrete path integral defined above and show its precise relation to the 
Ponzano-Regge model. For this purpose, we introduce the following parametrisation of SU(2) group elements 
in the fundamental two-dimensional representation (spin-^): 



G — COS01 + ismOu.a ~ e-\/l — jP- \ + ip.ff, 

where the matrices (Tfe, k = 1,2,3, are the (Hermitian) Pauli matrices normalised so that (cfc)^ = I for all 
/c's and Xxiauai) = l^ki- The parameter d e [0, 27r] is the class angle labeling the equivalence classes of 
group elements under conjugation (the rotation angle is 20), while u G labels the group elements within 
each equivalence class (it indicates the rotation axis). There is an obvious identification between the group 
elements g{6, u) and g{—0, —ii). We can thus restrict the range of the class angle to 6* € [0, tt]. The vector p 
is the projection of the group element G onto the Pauli matrices: 

p{G)^^MGa. (12) 

This vector has a bounded norm |p|2 < 1. The sign e = ± is the sign of cos9, that is e = + when 9 G [0, ^] 
and e = — when 9 G [f ,7r]. Finally, we also introduce the vector u = 9u which satisfies G = exp jiu.a). 
We can easily express the normalised Haar measure in term of these variables (see e.g. [tI I11||): 

dg . [\^Hd9 / ^ = E A / ^ = A / f^V d^^- (13) 
su(2) ^ io Js^ 4^ ti 27^' Jp<i J^^^ \ u ) 

Finally, we introduce the SU(2) characters Xj (G) for j e N/2 as the trace of G in the (spin-j) representation 
of dimension d^- = (2j -|- 1): 

, sin di9 

X,(G)-— ^-[/2,(cos0), 

where the functions C/„ are the Chebyshev polynomials of the second kind. 

We are now equipped to study the discrete path integral introduced above. Following 0, [HI) we get: 



(14) 



^ The maximal tree T touches every vertex v £ A but contains no loops; likewise for T* 
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Writing = p{Ge), we can integrate over the variables and obtain the product of d^^\pe) distributions. 
Now, the crucial point of the construction is the following identity: 

" g^B-plG) ^ 2^2^(3' (p(G)) = (5(G) + <5(-G) = d{G), (15) 
47r 

where 6{.) is the distribution localising the group element G on the identity in S0(3), i.e it does not 
distinguish I and —I as SU(2) group elements. Using the Peter- Weyl theorem, we can decompose S{.) into 
SU(2) representations: 

?(G) = 2^d,x,(G). (16) 

J en 

Then, expanding the 5{Ge) distribution in representations and writing explicitly the group elements Ge as 
the product of the g^-^ holonomies, one can recover the standard formula for the Ponzano-Regge model as a 
product of 6j-symbols attached to each tetrahedron of the triangulation A: 

{ieGN} e teA 

We insist that the sum is over integer spin representations je G N and thus this path integral formulation 
gives the SO (3) Ponzano-Regge model (as shown in Q). 

To recover the full SU(2) Ponzano-Regge model with a sum over both integer and half-integer representa- 
tions je € N/2, we introduce a new trick. The goal is to kill the S{—G) term in the 6{G) distribution. This 
can be achieved by slightly modifying the amplitude in the path integral: 

The new factor vanishes at G = — I and its value is simply 1 at G = I. Actually, we could use any other 
factor with the same property and the present one is the simplest such choice. Then, using the fact that 
tr(G) — Xi/2(G) by definition, one can check that: ^ 

1 [l + itrG] m = E ^^■^^(^) = -^(G)- (18) 

J6N/2 



At the end of the day, this final path integral with the extra factors reproduces exactly the Ponzano-Regge 
amplitudes for 3d quantum gravity including the sum over all odd and even representations of SU(2). 

B. Generalising the path integral 

Up to now, we have worked with an action S — -^i tr(i?eGe) = Be ■ p{Ge)- We could however use 
any other vector 0(G) instead of p{G) as long as the distribution S^^^cj)) coincides with 6{G) (or S{G)). We 



A standard formula in the recoupling theory of SU(2) representations gives for all j G N*/2: 

XjXL = X,-4_i + X,_ 1 • 

2 •'^2 J 2 

This leads to the following formula for any sequence of coefficients {oj}jgi^: 
We also point out the fact that Xj{~G) = (— l)^-'Xi(G). 
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also require that (p has a nice behavior under SU(2) conjugation: if we act on the group element G — > kGk ^ 
with k € SU(2), then the vector should simply rotate as a 3d vector (f> ^ kcj). 

Taking into account these two conditions, we introduce a whole family of possible actions: 

S = Y,Be-$iGe), with 0(G) = 0(0) Zi, 

e 

where 6 G [0, tt] t-^ (j){6) e M is a continuous map which vanishes only at = and possibly also at = tt. 

A subclass of actions is given by (continuous) positive functions of the vector p: = f{p)p with f{p)>0 
is a rescaling of p depending on its norm. The advantage of this subclass of actions is that p{G) is a nicely- 
behaved polynomial of the group element G. The particularity of this choice is that it is symmetric under 
^ TT — and will never distinguish I from —I. Thus it will naturally lead to the 5(G) distribution unless 
we use the extra factor which we discussed above. 

Other possibilities shall not be expressed solely in term of p and will in general also depend on the sign e (i.e 
whether is smaller or larger than 7r/2). A natural example is simply (f) ^ u ~ 6 u. Unfortunately it is not a 
polynomial of the group element G. There are other interesting choices such as (j) — — itr(G(7)/tr(G') ~ tanflw 
but we will not discuss them further.^ 

Since all these choices define the same path integral and lead to the same partition function (up to the 
5{.) vs (5(G) subtlety), the difference between these different formulations resides in the Lagrange multiplier 
i.e the _B-field. Indeed, the Lagrange multipliers differ by a G-dependent rescaling which will affect the value 
of observables depending on the B variables. More precisely, we write the generalised path integral using an 
arbitrary vector: 

where we replace the i?e variables by variables to underline that they do not represent the same field. In 
the case of a simple rescaling of the p vector, (p ~ f{p)p, we perform explicitly the change of variable and 
define Be = f{p)Xe. This gives: 

If we now insert in the path integral Z'-'^^ a X^ for some edge e, it will thus correspond to the insertion of a 
Bl/ f{pe)^ term which obviously depends on the group elements ge'- We will compute this effect explicitly 
in the next section when studying the length spectrum. 

Let us point out that each choice of momentum vector 0(G) defines a different Fourier transform between 
SU(2) and 'M? and induces a different covariant differential calculus on the non-commutative space jl7| . 

The issue behind this freedom in defining the discrete path integral can be seen from two perspectives. 
Either we can see it as the problem of identifying which discrete B^ variables are truly the discretisation of 
the triad field B in the continuum. Or we can consider it as an ambiguity in the definition of the measure 
for the discrete path integral. We do not address this issue in the present work. 



We can introduce the following vector for a group element G G SU(2) 

^ tr{Gcr) 



itxG 



tanSfi. (19) 



Contrary to the vectors p or n, the range of w is not constrained and sweeps the whole Mr . The normalised Haar measure 
reads in term of w: 

/ '^9= — TT—^- (20) 

This measure is actually of the type d^p^ /(I +p'^ /k^)" usually considered when deriving generalized uncertainty principle 
(GUPs) taking into account a minimal length in phenomenological models for quantum gravity. 
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III. B-INSERTIONS IN THE PATH INTEGRAL 



We would like to discuss the insertion in the discretised BF path integral of amplitudes depending on the 
B-field. Physically, this corresponds to either computing the average value of _B-observables or to taking 
into account a potential depending on the S-field. On the one hand, S-observables represent geometric 
quantities since the triad B defines the metric. On the other hand, introducing a _B-potential allows on 
to spinfoam quantise a whole class of "constrained topological theories" including Yang-Mills theories and 
gravity p^ . 

Here, we focus on (ultra-)local i3-insertions which depends on the i3-variables through the norm of the 
Be variables without any coupling between the various edges of the triangulation. The simplest case is the 
insertion of B^ along some edge e. Mathematically, it means that B^ is not anymore simply the Lagrange 
multiplier enforcing the triviality of the holonomy around that edge. Physically, this calculation gives the 
"length spectrum" for the spinfoam model. 

A. B^-insertion and the length spectrum 

Let us focus on a single edge e of the triangulation and insert B^ in the path integral along this edge. 
Instead of the original amplitude exp{iBe ■ Pe), we now have the weight B^ exp{iBe ■ Pe)- Integrating over 
Be does not give anymore the 5{Ge) distribution. More precisely, we compute: 



/ 



d'B^.B.,,G) ^ 2j2d,xAG), J ^^B^e^S.,(G) ^ 2^[4j(j + l)-3]d,x,(G). (21) 



47r 



We extract the length spectrum from this formula: 



3 /. 1^ ^ 



i ■ = J{J + 1) - i = + 2 1 ^ (22) 



First, we point out that Xq = —4/3 < for the trivial representation j = which defines a negative "vacuum 
area" and that — would vanish for j = 1/2. Second, this length spectrum actually holds only for 
integer representations j E N. To obtain half-integer representations, we need to work with SU(2). To this 
purpose, we insert the factor (1 + itrG)/2. This gives a slightly different and surprising formula: 

^ / ^ (1 + ^trc) e^^-f'^^^ = 5][4j(j + 1) - 3]d,xAG) + ^ 4j(j -I- l)d,x,(G), (23) 

where the —3 shift of the length spectrum appears only for integer representations. Notice that despite this 
shift, the length spectrum still increases with the spin j without anomaly. 

It is always possible to insert another factor F(G) depending on the holonomy, thus changing the measure 
on G, in order to shift arbitrarily the length spectrum There actually exists an infinite number of measure 
factor reproducing the same length spectrum. 

In order to derive these formula, we have to compute the integral / (PBdG Xj(G)B^ exp{iB ■ p). There 
are various ways to deal with it. We can first perform the integral over the group SU(2) and we obtain^: 

dGxAG)e'^-'~ = ^M\B[h (24) 



For instance, we consider an extra factor such as a -|- (1 — a)xi{G)/3 which is equal to 1 at G = ±1 and does not affect 
the partition function. Nevertheless, it changes the length spectrum for J S N leading to: = + 1) — {1 + 8a) /3. For 
a = —1/8, we use the measure factor {3xi(G) — l)/8 and obtain = + 1). On the other hand, we can use the measure 
factor (xi(G) — l)/2 with a = —1/2 and obtain simply Lj = d^. 
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where the J's are the Bessel functions. This gives the probability distribution of the -B-field and we get: 



/ 



d-'BdGxj{G)B^e 



2AB-P 



-I 



B^dB 



2Jd,{B) 
B ' 



Unfortunately, this integral does not converge. Thus, in order to compute this integral, we use the standard 
QFT method and introduce a source term exp(— ii? • q). B-insertions will be represented by differential 
operators of the type d/dq and we will send the source g to to obtain the final value. This gives for integer 
representations j G N: 



1 f d^B 
47r 



dGxj{G)B'e 



2AB-P 



1 



d d 



dqk dqk J i/l - q' 



(25) 



g=0 



where the 1 / -^/l — q'^ factor comes from the Haar measure and the group element g is given in term of 
the vector source q as g = — + iq.a. This is a straightforward Laplacian calculation in spherical 
coordinates 



A 



4 + 4rf,. = -dj{dj-2){dj+2). 



Keeping in mind that q = sin0, we compare this calculation to the standard Casimir formula when we act 
with the Laplacian on SU(2): 



" sin 



= -H-K- - l){dj + 1) = + 1). 



0=0 



B. Comparison with other path integrals 

It is natural to wonder how much does the length spectrum depend on the initial choice of path integral. 
The effect is comparable to ordering ambiguities in the quantisation process. It will not affect the leading 
order in j{j + 1) but will create a lower order shift compared to the previous calculation. This entirely comes 
from the difference in the Haar measure factor in the various parametrisation of the SU(2) group elements. 



^ We use the following spherical integral: 

-^2.-. - sin|X| 



Js2 47r 



/S2 47r \X\ 

We can also notice that the function to differentiate is easily expressed in term of c s cos = \/\ — cp". 

x(g) ^ Vii{<^) 

in term of the Chebyshev polynomial (of the second kind). The evaluation of the Laplacian Ag at </ = actually collapses to 
a first order derivative in c: 

N/(c)l,=o = -3/'(l), 

when /'(I) is finite. Differentiating the Chebyshev polynomials can be easily done many different ways. An elegant method 
is to use their generating functional: 
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For instance, let us consider the path integral defined with the vector u — 9 u and compute 
J (fiBdG B^Xj{C!) exp{iB ■ u). Then the corresponding length spectrum is given by the following Lapla- 
cian evaluation: 



d d sin^ u 

Since the norm of the vector is simply u = 9, this is equal to: 



(26) 



2 , 2 \ sur 



The difference between these results and previous length calculations [13, in the Ponzano-Regge model 
resides in the Haar measure factors. 

This difference in the path integral measure affects the value of all observables depending on the Be 
variables. There does not seem to be a unique choice of path integral and thus a unique "length spectrum" 
in this formalism. We insist that this issue is not resolved by requiring that the discrete path integral defines 
a topological state sum but it is an issue of the continuum limit: which choice of discrete Bg variables 
represent faithfully the continuous triad field B (from which the metric is reconstructed). 



C. Polynomials B-insertions and further 

We can follow the same method as above to account for any polynomial insertion of the _Bg" type. Using 
the vector source and inserting more derivatives, we obtain: 

-j—dgB\,[g)e^'^^^ ^ (2j + 3)(2j + l)(2j - 1), 

- y —dgB\,{g) e^*-'^s ^ (2j + 5)(2j + 3)(2j + l)(2j - l)(2j - 3), . . . (27) 

and so on. This sequence as n increases is actually related to the expansion of (modified) Bessel functions. 
Indeed, let us now insert in the path integral a Gaussian weight exp(— AB^) i.e a quadratic potential in B. 
Such an insertion allows the spinfoam quantisation of 2d Yang-Mills theory. We first perform the integration 
over Be'. 

Jd'Be-^^\^^-^^ = 

Then, without even having to introduce a source, we compute (see [l8l | for similar calculations on SU(2) 
Gaussian states): 

/^ci,e--\,(,)e^-^ = ^ f^^xMe'^ = ^e"^ L V ( ^) 
2j An ^ Ajvyy 8(^A)ii ^ 4Ai ['\8>^J ' \8^J. 

The previous polynomial insertions in _B^" can then be extracted from an expansion of this formula in A ~ 
(which corresponds to the asymptotic at infinity of the Bessel functions). 



D. Feynman propagator insertion 

In order to take into account particles propagating on the spinfoam, we can also insert Feynman propaga- 
tors along the edges of the triangulated space-time manifold following the framework introduced in [111 . Il2l | . 
This allows one to add particles to the Ponzano-Regge model directly at the path integral level. 
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We restrict ourselves to spinless particles with mass m — k sin (j) given in terms of the Planck mass k and 
an angle (j) S [0, Let us insert along a given edge e the (undeformed) Riemannian Feynman propagator, 
then the piece of the path integral concerning this edge reads: 

^ ^—i{sin (f}—ie)\B\ 

d^BeF„,iB,)e'''^ P', with F„,{B) = — , e ^ 0+. (28) 

At:\B\ 

The integration over B,, is straightforward and gives the standard momentum representation of the Feynman 
propagator as expected: 

d?B,Fm{B,)e'^^-^^ = ^2——. (29) 

— sm + le 



Following 12], this hmction on the SU(2) group is easily expandable in representations 
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KJG) = ^ = ^-V e-*'*^('^-")x7(G). (30) 

This defines the corresponding Feynman propagator insertion in the spin foam amplitudes. Moreover, it is 
given in the usual format, in terms of SU(2) representations. 

Such Feynman propagator insertions allow one to couple particles (and fields) to the spinfoam path integral 
and to study the coupled and effective dynamics of both matter and gravitational sectors [111, [l^ • 

The insertion of an on-shell particle corresponds to the imaginary part of this propagator: 



Atari c 



where 6{-) is the SU(2) distribution identifying the class angle of G to as a S0(3) group element, that is, 
without distinguishing the angles (j) and (tt — 0). Also, the imaginary part of the initial _B-insertion is: 

[F„,(|i?|)-f™(|i?|)l = ^ C0S| f ^2^^.B.si„,., 

L '^J 87r|B|tan<?!) 327r2 Jg2 



4tan( 



Written in this last form, the insertion of the massive on-shell particle is done directly at the level of the 
discrete action (as it was first shown in [3], up to the cos0 factor which was discussed in [Til. [20l|). 



This formula can be checked directly since it only involves geometric series. One can otherwise compute the coefficients of 
the decomposition into representations by evaluating the projection of the Feynman propagator onto the SU(2) characters; 



f 1 2 r 

J \p{g)\'' - sm^ (p + te n Jo 



(le - 



This integral is evaluated using the residue formula for j > 1 and e — > 0"*". It vanished for half- integer spins j G -I- ^ and 
is equal to 2e~"^j'^/ cos<f> for integer spins j G N. More precisely, we write the previous integral as a contour integral in the 
complex plane: 



dz 2"-l 



1 /■^'^ sinee'^j" 1 r 

— / 5 5 = . 

ii" Jo sin 9 — sin (p + ie tt Jc 



Ic (^2 - sin^ (ji + ieY 

where X = sin0 = {z — z~^)/2i and C is the unit circle in C. The denominator can be expanded and we find poles at 
±exp(iti</)): 

1 _ (2jz)2 

(X2 - sin^ (p + it) ~ (z- e^'i' - (1 + i)e){z + e-^'i' - (1 - i)e){z - e-^'i> + (I - i)e){z + e^'i' + (1 + i)e) ' 
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A last remark is about the SU(2) Feynman propagator instead of the S0(3) Feynman propagator used up 
to now. We define it, at the representation level, as the previous sum (pO|) extended to the whole of N/2: 



cos(j) k^{G) 



2 ^ g-id.j{4>-ie) 
jGN/2 



1 



2 sin^ 



sin^ ^ + ie 



cos 6 — cos (p — ie 



(31) 



We notice that summing over all half- integer spins cuts the angle by half and kills the SO (3) identification 
between the angles (fi and (tt — cj)). We can not implement this SU(2) Feynman propagator by a simple 
_B-insertion. Indeed, any simple _B-insertion in the original path integral will naturally give solely SO (3) 
structures since it only depends on the vector p{G) . We need an insertion depending non-trivially on both 
Be and the holonomy Ge'- 



k^AGe)^ / d''Befm{Be,Ge)e 



(32) 



with 



fmiB,G) EE 



2 cos ( 



-XiiG) 



-?(sin (p—ie)\B\ 

47r|S| 



(33) 



IV. B-OBSERVABLES AND GAUGE FIXING 



Having discussed how the insertion of B-dependent functionals affects the discretised path integral, we 
aim to evaluate the partition function with such _B-inscrtions and to obtain the value of such _B-observables: 



{0)^^J VloVBO{B,uj) 



where the normalisation is given by A/" = Zm. Both numerator and denominator need to be gauge-fixed. 

In the discrete setting, we choose a fixed edge e of the triangulated manifold A and we seek to evaluate 
the average value of the length (squared) : 

2a - n / <^'Be n / dge' e'-^lB.,,,.]^ (34) 
iBi) = ^ n / d'Be n / dffe. Ble^^^^-^^'\ (35) 

with the action S — ' P{Ge)- 

First, this requires proper gauge fixing |7i]. We choose a maximal tree T g A and we fix = for 
all edges e € T. As we will show in an example later, we actually require that e does not belong to T . 
Without such gauge fixing, we would always get a trivial (divergent) result. For this case, the Fadeev-Popov 
determinant is trivial, but as we show in Appendix [Bl for more complicated settings, this need not hold. 

Second, we introduce a background "state" ■(/'(Gg), normalised such that ijjil) = 1. Inserting such an 
observable in the path integral does not affect the partition function, i.e computing the value of (■(/'(Gg)) 
gives 1, since the path integral contains a (5(Gg). On the other hand, such a state allows one to locate the 
edge e: the state ip{Ge) excites different modes of the holonomy around the edge e and will modify the value 
of Bg. Thus we propose to compute: 

(B|V(Gg)> ^ ^ n / n / "^de- B|V(Gg)e^'5[^-'^=*], (36) 

for the choice ip{Ge) — Xj{Ge)/dj. The normalisation dj is to ensure that = 1. The choice j = gives 
the vacuum value (-B|). 
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This setting ensures that computing Xji^s)) leads to the same calculation as in the previous section. 
Indeed we find: 

},IBixAG,)) = }J ^dgB'x,l.9)e''-" (37) 

In particular, the vacuum value of Bg is surprisingly negative, = —3. This can be compared to a 

non-zero vacuum energy. What matters for the spectrum is the length variation: 

^ l-{Blx,iG,)) {BD = 4j(.? + 1). (38) 

Thus, we recover the standard length (squared) spectrum given by the SU(2) Casimir operator. 



A. The tetrahedral triangulation 



As an example to illustrate the gauge fixing procedure, we consider the triangulation of the 3-sphere 
with two tetrahedra. This triangulation contains four vertices, the six edges connecting them, the four 
triangles being these edges and the two 3-cells (inside and outside tetrahedra). The dual triangulation 
(or spinfoam) has two dual vertices (representing the inside and outside tetrahedra) and four dual edges 
connecting these two dual vertices and going through the four triangles. Finally, there are six plaquettes (or 
dual faces) transverse to the six edges of the triangulation; see FIG. [T] 




FIG. 1: The triangulation of S'^ (dashed lines) and its dual (full lines). The triangles are labeled by single numbers, 
while the edges are coloured by pairs. One can see immediately that the dual edges are in 1-1 correspondence with 
the triangles. Furthermore, each dual face (plaquette) consists of two dual edges. 



We will number the four triangles (or dual edges) by fc = 1, . . . , 4. We denote an edge by the (symmetric) 
couple of triangles that share it [kl). Then, the discretised partition function reads: 

Z = j \{d?BM /"n'^5fceA^'=<'*'-^-^"=sr\ (39) 

k<l •' k 

Let us consider the observable (-B12) defined as the quotient of the Bjj'inscrtion normalised by Z. In the case 
where we do not gauge fix the partition function, both numerator and denominator contain divergent ^(I) 
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terms. Although these "obvious" divergences cancel each other, we are left, nevertheless, with a meaningless 
(and divergent) correlation: 

We now gauge fix the discrete path integral by setting all Be variables to along a maximal tree in the 
triangulation. As maximal tree, we choose T = {(23), (34), (41)}. We then recover the result given above: 

/rf3B/[rf5]x4etrS,.3.-^(5i53-i)5(525ji) ' 

where we have not imposed the 6{-) constraints around the edges of the maximal tree. After integrating over 
the superfluous group elements, we are left with: 

~ Jd^BjdGe^<BG) ~ ^^^^ 

which we have evaluated previously. Computing {Bf2Xji9i92^)) docs not change the gauge fixing procedure. 

The last subtlety is the case where we include the considered edge (12) in the maximal tree T. Then, the 
gauge fixing is not enough and we once again obtain a meaningless result. 



B. The ©-graph and further S-observables 

In this final part, we generalise the previous Bg calculation to a coupled observable tr{BaBh) probing the 
correlation between two edges, a and b, of a triangle. As a source for the length, we use the following 
9-state to excite the holonomies around the two edges a, b: 

^^'^"'^'^ ^ dil'd'j^'-^^''^^''^^''^^'^^''^''^- ^^^^ 

Wc have normalised this state in order to ensure that '0(11, 1) = 1, so that its insertion in the partition function 
does not modify its evaluation. The representations ja^jb respectively excite the holonomies around the edges 
a and b while the representation J couples the two edges. Following the calculations of the previous section 
introducing vector sources qa,b for both edges, we find: 

d d 1 
{Ba-Bb) = - — — —— 'ip{g{qa),g{qb)) 
9qadqb^{l-^){l-^) 

This can be computed straightforwardly as well as the average values {B^) and {B^): 

(Ba-Bt.) = 4-1, (43) 
{BD = (4-1) + (4-1) -3, 
(4^) = (4-l) + (4_i)_3. 



Were we to consider two arbitrary edges of the triangulation, we would have to rotate the B-fields into the same coordinate 
system, using the holomony variables, in order to arrive at a gauge invariant observable. 
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When J = 0, then dj = 1 and we recover the previous resuh for the length of a single edge. We check 

that the scalar product {Ba ■ Bt) is correctly smaller than the product of the norms {B'^) (_B^) (as soon as 
ja > 1 or jt > 1). 



V. CONCLUSION 



In this paper, we presented a short note on geometric observables in 3d discrete quantum gravity. In 
particular, we concentrated on the accurate evaluation of polynomials in the _B-field, within the path integral 
approach. 

The discretised version of the BF theory path integral has the interesting property that is kills half the 
degrees of freedom of the SU(2) gauge theory (the representations j G N + i). It maps it to an S0(3) gauge 
theory. Fortunately, one can re- introduce them by inserting an appropriate observable (|17p. We noted that 
further ambiguity enters in the choice of Lagrange multipliers (J5 fields). 

By far our main result was to compute the expectation value of gauge invariant polynomials of the B field. 
We focused our attention on the length operator, and introduced the necessary gauge fixing to measure the 
observable. The gauge fixing occurred in two stages. The first was to introduce a normalised "state" which 
served to pick out the edge of interest and peak it on a certain representation. The second was the familiar 
gauge fixing of the isu(2) symmetry using maximal trees. We found that the expectation value of the length 
squared was actually negative for the vacuum state {j = 0) while it was greater than or equal to zero for 
all higher representations. A rather relevant fact is that the choice of discretisation of the B field affects 
the expectation value of the length variables. Thus, while the ambiguity has no net effect on the partition 
function, it has a very real consequence when observables are computed. 

We undertook the generalisation of this process to more complicated observables, such as higher order 
polynomials in the length, and other gauge invariant quantities, such as the dihedral angles between two 
edges. In particular, we could show that the edge vectors of a triangle satisfied the Cauchy-Schwartz relation. 

We elaborated on the topic of inserting an observable corresponding to the "Feynman propagator" of a 
point particle. Although, this concept is more an issue in the Lorentzian theory [l^l, observables of an 
analogous functional form can be introduced in the Riemannian regime. We found that we can define such 
propagators for both the S0(3) and SU(2) theories. Indeed, the "Hadamard" function, describing on-shell 
propagation, could be defined in the usual fashion: as the imaginary part of the Feynman propagator. 

Finally, our work here has revealed an avenue to deal with B observables in a more general context. The 
next step would be to compute more coupled i?-observables in order to study the discrete path integral for 
constrained BF theories with a i3-dependent potential such as Yang-Mills theory and gravity in more than 
three space-time dimensions. 

APPENDIX A: A FEW USEFUL FORMULAS 

A first interesting formula is the representation of the wave exp(ji? ■ p) (see e.g. [2l|): 

eiMB.) ^ ^ e-^d.^-^^^Ue'-'-'g), (A.l) 
where we have decomposed the vector B — \B \ B into its norm and direction. 



The simplest method is to use the leading order behavior of the character 

X^big}) ~ rfj l + cJ■ 



.2(l-'^?) 



15 



Secondly, we give the inverse formula allowing one to derive the characters from the wave and Bessel 
insertions: 

SB^-'-'J^ ,B sin(B|p1)J.,(B) . ^,,(«). (A.2) 

with IpI = sin^? > 0. The pre-factor is due to the Haar measure and is equal to — ~ ^Xi/2{g)- 
This defines the Bessel i?-observable that one has to insert in the path integral in order to induce the 
character Xj{G) on a plaquette instead of the basic 5{G). This can be used when we want to "gauge-fix" 
the representation je living on a given edge e to a fixed value. 



APPENDIX B: THE FADEEV-POPOV DETERMINANT FOR B-OBSERVABLES 

We are interested in computing the Fadeev-Popov determinant in more detail. For i?-observables, it 
emerges that it is non-trivial in many cases, although it was shown in 15] that the Fadeev-Popov determinant 
for pure gravity is trivial. 

Let us consider the triangulation A and the graph F, upon whose edges one inserts arbitrary _B-observables. 
F need not be connected but may consist of several components. The important point is that the B- 
observables are SU(2) invariant but not translation invariant. Thus, on a component Fj C F with \v\i 
vertices, the translation symmetry is broken at — 1 of them. This means that when one wishes to pick a 
maximal tree T C A, it should only intersect each component F^ at one vertex exactly. Furthermore, when 
dealing with a tree, one can always pick a vertex to be the root, and orient the edges so that they point away 
from this root vertex. In particular, it should not contain an edge of F. We shall neglect the gauge- fixing of 
the SU(2) symmetry as the observables are invariant and Fadeev-Popov determinant is trivial [l5|. 

Now that we have a maximal tree, we utilise the translation symmetry to set = Q for every e € T. We 
recall that the discrete version of the translation symmetry (fTO)) is: 



where we define the functions [/" and il" later on. 

The inverse of the Fadeev-Popov determinant for such a symmetry is: 

DfP= [U 6{Ul^'Utie) - K^'Kcl^tie)] - f/|^'^V.(e) + m^'K cj^^ie)]) ■ (B.l) 

•' eGT 

The Jacobian for the change of variables: 

0e ^ 0e = C/*^^Vt(e) " <Ptie)] " f/|^'^0s(e) + C^^ie)] (B.2) 

for every edge is: 

.({.,).{4))^ri „.,.,,(„.,.,;^^„,,.,.y (B.3) 

where fl — ii^ ■ a. Thus, the Fadeev-Popov determinant is: 

Dfp = n If^e^'^l {\Ul^'^f+m'}'^\') . (B.4) 

eGT 

Generically, the terms and fig are functions of the holonomies associated to the edges incident at v apart 
from the edge e itself. They have the property that — 1 and fig = when the curvature of all the other 
incident edges vanishes. 

This fact is enough to deal with the simplest case of a graph F with only one connected component. As 
mentioned above, the tree only hits F at one vertex, and we shall call this vertex the root Wroot- For the 
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moment, however, let us focus on a vertex at which the tree ends, that is, where only one edge of the tree is 
incident. Then, the quantum amplitude ensures that the curvature around each of the edges e ^ T incident 
at the vertex v vanish. Thus, the Fadeev-Popov factor for that edge is trivial. Furthermore, the Bianchi 
identity ensures that the curvature associated edge e G T incident at v also vanishes. By applying that 
procedure to the outer edges and working our way towards the root of the tree, we can show that the whole 
Fadeev-Popov factor is trivial. 

For the case of a graph F with two components, we call the vertices where the tree T intersects Fi, F2 as 
■yroot and utcrm respectively. In particular, there is a unique path p C T joining Wioot and wtorm- When we 
attempt to apply the same procedure as above, we find that the Fadeev-Popov factors are trivial for edges 
e £ T/p. But for edges e S p, we find that the curvature associated to the two edges e G p incident at a 
vertex are non- vanishing. Thus, the Fadeev-Popov factor for these edges factor non-zero. 

We shall calculate [/" and fig explicitly in the case v £ p. We know that the translation symmetry is 
satisfied thanks to the Bianchi identity: 



n 



9e = 



(B.5) 



where Qe is the holonomy associated to the edge e. But the curvature vanishes for all but the two edges 
e G p, so the identity reduces to: 



9192 = (Mil + Pl)(u2l + P2) = (^UiU2 + itr(PiP2) j I + (^UiP2 + U2P1 + ^[Pi, P2] I = I. 

Then, the transformation of the relevant B variables is: 



(B.6) 



Bi Bi+U20+-[P2, 
B2 B2 + UiC^ - ^[Pi, 



Say that ei terminates at v. The Fadeev-Popov determinant for that edge is: 

ic/n(ic/rp+4|fi5:p) = Ki(Kp + ^iP2 



1 



cos O2 1 I cos^ ^2 + ^ sin^ 62 



(B.7) 
(B.8) 

(B.9) 



But the Bianchi identity imposes that 9i = 
finishes up as: 

Dfp = 



92 ='■ 0. Therefore, the Fadeev-Popov determinant for this tree 

iPl 



cosy cos 



1 . 2 

- sm 

4 



(B.IO) 



where we li3,ve assumed tticit the tree is cliosen such that i^tcim 

lies as an endpoint of r2. Should it lie in the 
interior, the Fadeev-Popov factor for the edge e G T incident at Wtcrm will be more complicated since more 
than two edges at Wtcrm will have non- vanishing curvature. 

For the more general case of graphs with multiple components, the Fadeev-Popov factors are trivial on all 
the edges apart from those on the unique path in T which join the vertices in T n F. Also, apart from the 
case where F has just one component, the Fadeev-Popov determinant, and hence the resulting amplitude, 
depends on the choice of tree. 
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